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SUMMARY 

An objective method utilizing least squares is presented for the 
determination of the optimum parametric constants for stress-rupture data 
The method is applied to both isostress and isothermal data for the pa- 
rameters proposed by Larson and Miller, Manson and Haferd, and by Dorn. 
Several examples are treated in detail, and it was found that the method 
gives good results. It is shown that the values of the constants for the 
parameter proposed by Manson and Haferd are not critical as long as T 

and log t a appear in the proper combination. In addition to optimiza- 
tion, the chief utility of the method lies in the fact that it gives the 
same results for a given set of data no matter who makes the analysis, 
which is not the case for the graphical methods presently employed. 


INTRODUCTION 

The problem of extrapolating and correlating high- temperature creep- 
and stress-rupture data is presently receiving considerable attention. 

In particular, the ability to extrapolate short-time rupture data can 
greatly reduce costly experimental programs and also reduce time delays 
in choosing a suitable material for a given application. 

The most widely used techniques at present for extrapolating stress- 
rupture data are the so-called parameter methods. These methods assume 
that by plotting the creep- rupture data for a given material in an ap- 
propriate parametric form a single master curve results which can then be 
used for interpolation and extrapolation purposes. The three best known 
parameter methods are those of Larson and Miller (ref. l), Manson and 
Haferd (linear parameter) (ref. 2), and Dorn (ref. 3). 

The Larson-Miller method assumes that a plot of log rupture time 
against the reciprocal of the absolute temperature at a given constant 
nominal stress is linear. Furthermore, it is assumed that all such 
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constant-nominal-stress lines intersect at a. common point (0,-C) (fig- 
1(a)). A plot , therefore, of the stress a against (log t + C)(T + 460) 
should produce a single master curve valid Tor all stresses and tempera- 
tures as shown in figure l(b). Originally, the constant C was taken 
equal to 20} however, it is now generally recognized that for best re- 
sults it will vary with the material. 


The linear parameter method assumes that a plot of log rupture time 
against temperature in degrees Fahrenheit at a given constant nominal 
stress is linear and that all such lines converge to a common point 
(T a , log t a ) as shown in figure 2(a). The constants T a and log t a 

are determined from the data for a given material in a given condition. 


T - T 


a 


should 


'a 


A plot, therefore, of the stress a against t _ log t 

produce a master curve, as shown in figure 2(b), valid for all stresses 
and temperatures . 


The Dorn method assumes that a plot of log rupture time against the 
reciprocal of the absolute temperature at s given constant nominal stress 
is linear and that all such straight lines are parallel with slope D as 

2.5 D 

shown in figure 3(a) . A plot, therefore, cf the stress against te 
should produce a single master curve valid for all stresses and tempera- 
tures (fig. 3(b)). 


It is seen from the foregoing discussion that, in order to make use 
of these parameter methods, certain materiel constants must be determined 
To do this, the general practice has been to plot creep-rupture data for 
a given material as shown in figure l(a), 2(a), or 3(a), depending on the 
selected parameter. The desired constants can then be obtained visually 
by appropriate extrapolation. Only when dsta are given in constant- 
nominal- stress form, can these plots be obtained directly. In the more 
usual case where the data are isothermal, several cross plots must first 
be made before figures such as l(a), 2(a), and 3(a) can be constructed. 


Since creep-rupture data generally have appreciable scatter, it is 
apparent that the results of such visual ci oss-plotting will depend upon 
the judgment of the individual analyzing the data. Furthermore, some 
experience in the field of material evaluation is necessary m constructing 
the various plots . 

This report presents an analytical method for determining the best 
values of the constants for the three parameters discussed. The method 
which is based on the standard least-squares procedure makes use of the 
original raw data and requires no judgment on the part of the analyzer. 
Moreover, since it is a least- squares procedure, it gives the statistically 
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most probable values of the desired constants. The method is first pre- 
sented for the case of constant-nominal- stress data and then for the more 
usual case of isothermal data. Examples are presented for several mate- 
rials for each of the two cases. 

It is not the object of this report to discuss the merits of the var- 
ious parameters used. This has been discussed at length in references 4 
and 5 where it is indicated that in general the linear parameter gives 
better agreement with experiment for extrapolated stress-rupture times 
than either the Larson-Miller or Dorn parameters. 


ANALYSIS 

Constant-Nominal-Stress Data 

Larson-Miller parameter . - Consider a set of constant- stress data as 
shown in figure 1(a). On the basis of the Larson-Miller parameter, it is 
assumed that a set of straight lines intersecting the log t axis at 
log t equal to -C can be fitted to these data. The equation of the 
straight line passing through the data for the first constant-stress line 
can be written as 


yd) a -c + b^xd) (la) 

and for the second isostress line, 

y( 2 ) a — C + b 2 xd) (lb) 

Thus, for any isostress line, say the j"^, 

y(j) a -C + bjT^j) ( lc ) 


where 


y = log t 

bj= slope of j^* 1 isostress line 

T = 

T + 460 


and the superscripts designate the particular constant- stress line under 
consideration. 





To find the best set of lines fitting the data and intersecting at 
the point -C, the sum of the squares of the deviations S of the actual 
data points from the lines (the residuals) is minimized. Thus, 



where p is the number of isostress lines and n^_, n 2-> and so forth are 
the number of data points for each line. 

In order to find the values of C anc the bj that will make S a 
minimum, S is differentiated with respect to C and the b^, and the 
resulting equations are set equal to zero. This results in 
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where 



and n is the total number of data points: 

n = + n 2 4* . . . + rip 

Solution of equations (3) for C and b i gives 

J 



By means of equations (5), the best value of the Larson-Miller pa- 
rameter C can be directly computed for a given set of constant- stress 
data. The best lines intersecting at -C can also be plotted by using 
the second of equations (5). 

Dorn parameter . - To determine the best slope D of a set of paral- 
lel constant-stress lines as shown in figure 3(a), a similar procedure 
as for the Larson-Miller parameter is used. The equation of any one of 
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the isostress lines can be written as follows: 

y( j) = a, + Dt( )) 

J J 

The expression to be minimized now becomes 


£ ^’T *£ [. 


(2l 

*1 " a 2 


”1 ^ 

4 2) ] * 



a - = minimum 

P 1 J 


Differentiating with respect to the a j and D and setting the resulting 
expression equal to zero give 

n .a • + A^D = ® j ) 


g £ 


B, = V C, 


where 


Bj_, and are as defined in equations (4) and 


B o - 2^1 


(8a) 
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The solution of equations (8) gives 



The Dorn parameter D can, therefore, be calculated directly from equa- 
tions (9) . 

Linear parameter . - A similar approach as in the preceding can be 
used for determining the constants T a and log t a (fig. 2(a)) for the 

linear parameter. The equation for any of the isostress lines is 

= y a " b j T a + b j T ^ (10) 

where 

Equation (lO), however, is nonlinear in the unknown constants because 
of the term bjT a . Minimizing the sum of the squares of the deviations 

would, therefore, lead to a set of nonlinear algebraic equations which 
would be very difficult to solve. To avoid this difficulty, two alternate 
approaches have been used. In the first approach, a value is assumed for 
T & . Equation (10) is then written as follows: 

y (j) = y a + b j JtC J ) - T 1 (lOa) 

Equation (lOa) has exactly the same form as equations (l) and the 
solution is, therefore, given by equations (5) with C replaced by -y 

and x replaced by T - T a . Once the best values of y a and bj are 
found for the assumed value of T a , the sum S is computed by equation 
(2) (again replacing C by -y a and x by T - T a ) . A new value of 
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T a is chosen, and the calculation is repeated giving a new value of S. 
The value of T a for which S is a minimun is the correct value. Since 

the results for the linear parameter are generally insensitive to the 
exact value of T a as long as the corresponding value of y a is used, 

the previous trial- and- error procedure need, in general, be carried out 
only a few times in order to obtain a satisfactory value for T a . 

An alternate simpler approximate method which does not involve trial 
and error can also be used. In this method, the nonlinear term b^T a is 

temporarily grouped with y a , and equation (lO) is written as follows: 

y(j) - d, + b.T^i) (11) 

J d 


where 


d j = y a - *j T a 

Equation (ll) is now linear in the unknown constants , dj and bj, and 

these can be found by least squares as before. Thus, the sum S to be 
minimized now becomes 


n- 


* d [4 x) - d i - h 1 

i=*l 




(2 ' - d 2 - b 2 T i 2) ] 2 + • ' ' + 


i=l 


- minimum 


- s - vi p) ] ! 

i=l 

Differentiating with respect to the dj anl bj leads to 


(12) 


n o d j + A j b J * D j 


A j d j + B j b j c j 


where A„ B., C,, and D, are as previously defined and 
J J V ** 


T . = T- 


( 13 ) 
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Solving equations (l3) gives 


(14) 


Since the best values for dj and b j have been determined from equa- 
tions (l4), the best values of T a and y a can now be found as follows: 


and it is desired to find the best values for y a and T a . The following 
sum is, therefore, minimized: 



which gives 


Therefore, 



(15) 


(16) 


( 17 ) 
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Thus, T a and y a are directly computed from equations (14) and (l7). 

It will be shown later that the values of T a and y a obtained using 

equations (l4) and (l7) differ very little from the values obtained by 
the trial- and- error method previously described. This procedure is more 
fully discussed in the appendix. 


Isothermal Dal a 

It has been the general practice to perform stress-rupture tests at 
constant temperature and to represent these- data as shown in figure 4. 
For this case a least- squares method to obtain the parametric constants, 
similar to the one previously described, can be used. To do this, the 
master curves shown in figures l(b), 2(b), and 3(b) will be presented in 
the following form: 

Larson-Miller parameter: 


(y 4 * C) (T 4 * 460) — aQ 4 * ci^x 4 * a^u.^ 4 * ... 4 * a^x^ 


Linear parameter: 


y - y a 

T - T a 


a r 


4 - a-, x 4 * aoX^ + 


• + ^ r 


(18) 


Dorn parameter: 


D 

y " T + 460 


a 0 + a~j_x 4 * 4 - 


• + ^ J 


vhere x S log cr, m is the degree of the polynomial assumed, and y has 
previously been defined. Note that the coefficients a-^ do not, of 

course, have the same values in the three equations (18) . 

In equations (18), the master curves are represented by polynomials 
in log a. In general, a parabola or cubic will describe the master 
curve with sufficient accuracy except for uhose materials whose master 
curve has a reversal in curvature near the tail end at low stresses. 
Generally, it is desirable to assume a parabola or a cubic first to deter- 
mine the constants as described herein. The master curve should then be 
plotted using these constants. If it appears that a reversal in curvature 
is present at low stresses, the calculation should be repeated omitting 
these data. Omission of these data should not affect the true value of 
the parametric constants since these are determinable from any segment 
of the master curve. It might be expected that increasing the degree of 
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the polynomial will improve the values of the constants. However, it will 
be shown that polynomials higher than the third degree are generally 
unnecessary. 

The least- squares method will now be applied to equations (18) . The 
quantity that will be minimized is the sum of the squares of the differ- 
ences between the logarithm of the actual rupture times and the theoreti- 
cal rupture times, since the rupture time is generally the critical 
variable . 

Larson-Miller parameter . - The first of equations (l8) is written 
as follows : 

y = -C + a Q T + a-jTx + a^Tx 2 + . . . + a^x™- (19) 

where t = l/(T + 460) . 

The sum of the squares of the residuals is 


S = 


2j {yi 


+ C 


a 0 T i " 


a l T i x i 


a 2 x i*l 


- ( 20 ) 


i=l 


Minimizing the sum of the squares of the residuals as previously done 
leads to the following set of equations for C and the a^: 


-nC + ®0 a 0 ^l a l 


Em 8 ® 


n 

E y i 

i=l 


-EqC + FQaQ + Fj^a-^ 
-EqC + Fj_aQ + F^q^q + * 

I / | 

i I i 

_E m C + F m a 0 + F m+l a l + 


• + Em 8 ™ 


J 0 




+ Em+l 8 ™ = G 1 

I I 

f I 

! | 

+ F 2m a m = G m 


( 21 ) 
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where 



The solution of equations (2l) gives the optimum value of C as well as 
the aj for an assumed degree polynomial for the master curve, 

A word of caution is needed with regard to the solution of equations 
(2l). These may in many cases be somewhat ill-conditioned, that is, a 
large number of significant figures may be lost during the process of 
solution. Thus, care must be exercised that enough significant figures 
be carried in the calculation to ensure meaningful answers. The situation 
is more aggravated the higher the degree of the polynomial assumed for 
the master curve. However, as will be shown, it generally would not be 
necessary to assume more than a cubic for nhe master curve. Note that 
the number of equations to be solved is equal to the degree of the poly- 
nomial assumed plus two . 

Linear parameter . - The second of equations (18) is written as 
follows : 

y = y a + (T - T )a n + (T - Tja.x + (T - T, )a Q x 2 + . . . + (T - T )a ^ 

^ cx \j a ± c. c a m 

(23) 

Equation (23) is of the same form as equation (l9) with -C replaced by 
y a and t replaced by T - T a . The solul ion, therefore, is given by 

equations (2l) and (22) provided T a is krown. A trial- and- error pro- 
cedure is thus followed as described for tire constant- stress data. Values 
are chosen for T a , and for each value equations (2l) and (22) are solved 

and S is computed from equation (20) (by replacing C by -y and 
t by T - T a ) . The value of T a for which S is a minimum is the 
correct value. 
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Dorn parameter . - The last of equations (18) is written as follows: 
y = Dr + a Q + a-j_x + a 2 x 2 + . . . + a^jX 111 (24) 

where t = l/(T + 460) . 

Minimizing the sum of the squares of the residuals leads to the following 
equations : 


naQ 4* EqD 4* Hj_a^ 4* ^2^2 


E 0 a 0 + F 0 D + E l a l + E 2 a 2 + 


H l a 0 

+ 

E X D 

+ 

H 2 a l 

+ 

% a 2 + ' 

H 2 a 0 

4 - 

e 2 d 

+ 

% a l 

4* 

H 4 a 2 + * 


I I \ \ 

II II 
//II 

E m a 0 + + H m+l a l + H m+2 a 2 + 


+ Vro = £ n 

i=l 

• + = G 0 
• + E m+l a Hl = I 1 y 


+ = I 2 

\ 

l 

\ 


+ H 2m a m “ 1 


m 




where Ej and Gq have previously been defined and 

n 

-E*i 


H 


i=l 


s E jA 


(25) 


(26) 


i=l ' j 

The solution of equations (25) gives the best values of D and the 
aj for a given degree polynomial for the master curve. 


EXAMPLES 

Constant-Nominal- Stress Data 

Consider a set of constant- stress data as given in table I and taken 
from a tabulation of the data in reference 4. It is seen from the table 
that 

n-j_ =9, r ±2 53 9, = 9, n^ = 8, n = 35 
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Letting X ± = — B — (eqs. (4)) gives 


A ± = 5. 1709x10" 3 

B i 

= 2.9758X10" 6 

C 1 

= 8. 0785x10“ 3 

= 13.748 

A 2 = 5. 299 7x10" 3 

B 2 

= 3. 1311x10“ 6 

C 2 

5.6297xl0“ 3 

D 2 = 8.8652 

A 3 = 5. 5102x10" 3 

B 3 

= 3.3793X10" 6 

C 3 

^ 4. 3022x10“ 3 

D 3 = 6.6343 

A 4 = 5.3898*10" 3 

B 4 

= 3.6358 xlO -6 

c 4 

= 5. 3599x10" 3 

D 4 = 7.6130 

and 

35 

T. 

y 4 = + D 2 + 

d 3 + 

I 4 =« 36.861 



i=l 

Then, from equations ( 5 ) 

C = 23.8 

b x = 44,146 b 2 =* 42,154 b 3 = 4C,151 b 4 = 36,819 

The best value of C in this case is therefore 23-8. 

By using equations (9), the Dorn parameter is calculated as 

D = 41,470 

To calculate the constants for the linear parameter by the trial - 
and-error method described, values must be assumed for T a - If it is 

not known at all what region T a might be ii, a rough plot could be made 
to determine a first guess for T a - Thus, assuming T a « 500 and, in 


this case, defining x = 

T - T a 

(eqs. (4)) g 

Lve 

A x = 7,049.99 

II 

i— 1 

PQ 

5.56625X10 6 

Cj = 10.2311X10 3 

A 2 = 6,694.99 

b 2 = 

5.06842xl0 6 

C 2 = 5.39997X10 3 

A 3 = 6,084.98 

b 3 = 

4.15518X10 6 

C 3 = 3.84065X10 3 

A 4 = 4,209.98 

11 

pq 

2. 23945 xl 0 6 

C 4 = 3.48307X10 3 
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and, as before, 

35 

y ± » 36.8608 

i*l 

Therefore, from equations (5), by letting y a take the place of -C, 
there is obtained 

y a - 11.4 

bp = -0.001266 b 2 = -0.001405 b 3 = -0.001583 b 4 = -0.001996 
Then, from equation (2), 

S = 0.6707 

Wow, by assuming a value of T a equal to 600, the computations are re- 
peated giving 

y a * 9.896 
S as 0.5819 

and, for T a = 700, 

y a = 8.306 
S » 0.5889 

These calculations indicate a minimum value for S at T = 600. If 

a 

greater accuracy is desired, a few more points can be taken in this 
vicinity. Thus, it turns out that the true minimum is approximately at 

T a => 650 

y a - 9-ios 

S = 0.5640 

However, as will be shown later, using values of T a = 600, y a =» 9.9 or 
T a = 700, y a = 8.3 will not affect; the stress-rupture results appreciably. 
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This trial-and-error procedure can he completely avoided if the 
approximate method described for the linear parameter is used. Thus, by 
defining t = T, equations (4) give 



11,550 

B 1 " 

14.866xl0 6 

C x = 17,105 

D x = 13.748 

11,195 

B 2 = 

14.013xl0 6 

C 2 » 9,832.6 

Dg = 8.8652 

10,585 

b 3 = 

12.490xl0 S 

C 3 = 7,157.8 

D 3 = 6.6343 

8,210 

b 4 = 

8.4494xl0 6 

C 4 = 7,289.6 

D 4 = 7.6130 


and from equations (l4), 


d x = 17.4 
d 2 = 17.9 
d 3 = 19.3 
d 4 a 23.4 

Substituting into equations (17) 


b x = -C .0124 
b 2 = -C .0136 
b 3 = -C.0158 
b 4 » -C.0219 
gives 


T = 643 
a 

y =9.3 
J a 


In this case the trial-and-error procedure is not really necessary as 
good results can be obtained by using equations (l4) and (17) to get T a 
and y a . 


Isothermal Datt 

Consider the isothermal data of figure 4 tabulated in table II. To 
calculate the best parameter C, let T = — + 4 ~ and assume a third- 

degree polynomial (m = 3) for the master cuive. Then computing the Ej, 
Fj, and Gj by equations (22) and substituting into equations (2l), with 
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m = 3, give the following five equations: 

-32C + 0.017013000 a Q + 0.067135551 & 1 + 0.26671741 a 2 + 

1.0666639 a 3 = 62.095000 

-0.017013000 C + 9.1227245X10" 6 a Q + 36. 213995X10" 6 a x + 

144. 71013x10" 6 a 2 + 582 .01515X10" 6 a 3 = 0.033195744 
-0.067135551 C + 36 . 213995x10" 6 a Q + 144. 71013X10" 6 a-L + 

582. 01515x10" 6 a 2 + 2355 .6764x10" 6 a 3 = 0.12753108 
-0.26671714 C + 144 . 71013x10" 6 a Q + 582. 01515X10" 6 a-[_ + 

2355 .6764x10” 6 a 2 + 9593 . 2956X10" 6 a 3 = 0.49259087 
-1.0666639 C + 582 .01515x10” 6 a Q + 2355 . 6764X10” 6 aj + 

9593. 2956x10" 6 a 2 + 0.039301429 a 3 = 1.9127579 

The solution of these equations gives 

C = 15.3 

a Q = 182,925 a-]_ = -92.231.9 a 2 = 20,404.8 a 3 = -1,697.03 

To calculate the linear parameter constants T a , y a , let = T^ - T a 
and replace C by -y a . For an assumed value of T , E^, F*, and G* 

are computed as before from equations (22). Thus, with m = 3 and 
T & = 0, equations (2l) become 

32y a + 46,000 a Q + 178,945.27 a ± + 700,948.43 a 2 + 2,764,570.3 a 3 = 62.095 
46,000 y a + 67.16xl0 6 a Q + 259 .11478xl0 6 a x + 1006 .6607X10 6 a 2 + 
3937.9381xl0 6 a 3 = 88,481 

178,945.27 y a + 259.11478x10 s a Q + 1006 . 6607xl0 6 a x + 3937.9381xl0 6 a 2 + 

15, 510.467x10 s a 3 = 335,601.61 
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700,948.43 y a + 1006 .6607xl0 6 a Q + 3937.93Elxl0 6 a ± + 15,510. 46 7 xlO 6 a g + 

61, 505 .824xl0 6 a 3 = 1,280,135.8 

2,764,570.3 y a + 3937 .9381X10 6 a Q + 15,510.467xl0 6 a 1 + 61, 505 .824xl0 6 a 2 + 

245,526.23 a 3 = 4,920,487.4 

The solution of these equations gives a value of y a = 15,0. The value 
of S obtained from equation (20) is 0.4751. Choosing a value of T a 

equal to 100 and repeating the previous calculations result in a value of 
y a = 14.1 and S = 0.4801. A minimum value of S = 0.4671 is obtained 

for T a =a -300 with a corresponding y a = 17.7. A plot of S against 
T a and the corresponding values of y a is shown in figure 5. It is seen 
that S does not change much with T a in the range of 0 to -500. It 
would, therefore, not make much difference which value of T a in this 
range is used as long as the corresponding value of y a is used with it. 
For other materials, the curve of S plotted against T a might have a 
sharper minimum, and the value of T a would be more critical. 

For the Dorn parameter, Ej and Gq are computed as for the Larson- 
Miller parameter, and Hj and Ij are computed using equations (26) 

^with ^ . Then, for m = 3, the following equations are 

obtained: 

32a 0 + 0.017013000 D + 125.50874 ip + 495.64303 a 2 + 

1970.5813 a 3 « 62.095000 

0.017013000 a Q + 0.0000091227245 D + 0.057135551 ap + 0.26671741 a 2 + 

1.0666639 a 3 « 0.033L95744 

125.50874 a Q + 0.067135551 D + 495.64503 ap + 1970.5813 a 2 + 

7886.7697 a 3 = 237.?6461 
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495.64303 a Q + 0.26671741 D + 1970.5813 a ± + 7886.7696 a 2 + 
31,770.400 a 3 = 911.60106 

1970.5813 a Q + 1.0666639 D + 7886.7697 a ± + 31,770.400 a 2 + 
128,793.83 a 3 = 3521.6275 
Solving these equations gives 


D = 32,900 


RESULTS AND DISCUSSION 
Determination of Master Curves 

In order to determine how well the various parameters are determined 
by the least- squares method presented, master curves for the linear pa- 
rameter were computed for the examples described. These master curves 
are plotted for the constant- stress data of the 17-22A(s) steel and the 
isothermal data of the 18-8 stainless in figures 6 and 7, respectively. 
These curves were then used to replot the constant- stress lines for the 
17-22A(s) and the isothermal curves for the 18-8 stainless. The results 
In figures 8 and 9 show that the agreement between experimental data and 
the computed lines is good. 


Effect of Polynomial Approximation 

As a further check on the assumption that the master curves can be 
represented by polynomials as given in equations (18) , the parameters for 
the constant- stress data of the 17-22A(s) were recomputed using a poly- 
nomial representation for the master curve. Thus, equations (2l) were 
solved for the linear parameters assuming polynomials of the third degree 
and also of the fifth degree for the master curves. The values obtained 
for both the third-degree and the fifth-degree polynomials were 23.8 for 
C, 650 for T a , and 9.1 for log t a . These values are seen to be the 

same as previously obtained using the constant- stress data directly to 
find the best intersection of the straight lines. 

As a further illustration, the data for 25-20 steel given in refer- 
ence 2 were considered. The parameters for this material were obtained 
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using second.-, third-, fourth-, and fifth-degree polynomials. The fol- 
lowing results are obtained: 


m = 

2, 

c = 

14.4, 

T a = 

-400, 

log t a = 

17.3 

m = 

3, 

c = 

14.2, 

?a = 

-200, 

iog t a = 

15.3 

m = 

4, 

c = 

14.4, 

T a = 

-200, 

log t a = 

15.3 

m = 

5, 

c « 

14.6, 

T a = 

-200, 

!°g t a = 

15.6 


It is seen that the same results were obtained for all the polynomials 
except for the polynomial of second degree <rhich gave slightly different 
values for the linear parameter constants. However, there is very little 
difference in the sum of the squares of the deviations in going from 
T a = -200 to T a = -400 so that even using a value of m = 2 would 

give good results in this case. 


Insensitivity of Results to Parameter Values 

As a further illustration of the relat Lve insensitivity of creep- 
rupture-data correlation to the precise values of the linear parameter 
constants (as long as T a and log t a appear in the proper combination), 

the data for Nimonic 80A as given in table III were analyzed. These data 
were taken from reference 5 where the parametric constants are given as 
16.9 for C, 660° F for T a , and 9.65 for log t a . In a private communi- 
cation to the authors of the present paper other investigators questioned 
the values given in reference 5 stating than their analysis using the same 
data gave values of T & = 100° F, log t & = 16 . 

An analysis was therefore made of these data using the least- squares 
method presented herein. The results are shown in figures 10 and 11. 
Figure 10 shows that the best values for the parameters are T = 400° F, 

cl 7 

log t a = 12.2. However, because the curve js flat in the region of its 

minimum (it is drawn here to a very expanded scale in order to show the 
precise minimum), other combinations of constants show sums of deviations 
not much higher than those at the minimum. Thus, the sum of the squares 
of the deviations S for the minimum point is in the neighborhood of 

O 

5.5x10 , the value for the constants of reference 5 is approximately 

- 2 

6x10 , and for the constants of the private communication it is approxi- 

mately 5.7xl0“ 2 . Thus, all three combinations are for all practical pur- 
poses equally good. This is further illustrated in figure 11 where the 
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reconstructed isothermals for each of the three combinations are compared 
to the experimental data. All the computed values lie on the same solid 
lines. It should be noted that all three combinations lie on the 
straight-line plot of log t a against T a of figure 10. 

This illustration shows graphically how different Investigators 
analyzing the same data can arrive at different values of the constants 
in the linear parameter unless an objective method such as the least- 
squares method presented herein is used. It also shows, however, that 
the degree of correlation of the data is rather insensitive to the pre- 
cise values of the parameters as long as the proper value of log t a is 

used with T a . This insensitivity of the correlation to the precise values 

of the constants is due to the fact that the intersection point of the 
constant- stress lines is generally remote from the actual data points. 
Therefore, moving the intersection point along an average line through 
all the data would not appreciably change the individual lines. It has 
been the experience of the authors that this is true for most materials. 


CONCLUSIONS 

An objective least-squares method has been presented for determining 
the optimum values of creep-rupture parametric constants for the Larson- 
Miller, linear, and Dorn parameters. From the examples shown it is con- 
cluded that the results obtained are insensitive to the degree of poly- 
nomial assumed for the master curve and that for the linear parameter the 
actual values of T a and log t a are not critical as long as they appear 

in the proper combination. Furthermore, the method permits a person with 
no experience in the field of materials to obtain the correct values of 
the parametric constants from tabular data. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, December 12, 1958 
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APPENDIX - APPROXIMATE METHOD FOR CONSTANT-STRESS DATA 

Consider a set of constant-stress data as shown in figure 8. The 
trial- and- error method for obtaining the best intersection point for these 
lines utilizes the least-squares procedure by minimizing simultaneously 
the squares of the deviation of all the data points from the lines. By 
this method each point is given the same weight. 

The approximate method described in the body of this report first 
treats each set of constant- stress data separately and finds the best 
fitting straight line for the set. Now, if the set of straight lines 
thus determined is to intersect at a common point, a plot of slope against 
intercept for these lines should be a straight line as can be seen from 
equations (ll) where bj is the slope of the line and dj is the 

intercept. The so-called "best average 1 ' intersection point is therefore 

found by fitting the best straight line to the plot of d^ against b^. 

J J 

It is to be noted that by this approxinate procedure each constant- 
stress line is given the same weight without regard to the number of data 
points associated with that stress. If the data are such that one or more 
of the lines is ill defined then it may be desirable to employ some 
weighting procedure when using the approximate method. 
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TABLE I 


STRESS-RUPTURE DATA FOR 


17-22A(S) (REF. 4) 


a 

T 

t 

a 

T 

t 

o 

o 

o 

•\ 

o 

1 — l 

1370 

2.8 

40,000 

1285 

0.075 


1370 

3.7 


1260 

.37 


1350 

4.5 


1210 

1.35 


1315 

12.5 


1210 

1.90 


1270 

48.5 


1175 

6.60 


1270 

51.3 


1150 

13.6 


1235 

129.8 


1120 

39.5 


1210 

228.7 


1100 

83.0 


1160 

1301 


1075 

205.7 

20, 000 

1400 

0.1 

O 

o 

o 

o 

00 

1140 

0.033 


1375 

.12 


1070 

1.5 


1320 

1.0 


1045 

2.5 


1270 

3.9 


1030 

5.7 


1230 

13.3 


1000 

15 


1190 

48.0 


985 

82 


1170 

102.7 


970 

109 . 2 


1140 

242.1 


970 

433 


1100 

987 
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TABLE II. - 18-8 STAINLESS STEEL (REF. 2) 


Tempera- 
ture, T 

log t 

a 

log a 

Tempera- 
ture, T 

log t 

a 

log a 

1200 

0 

38.0*10'-’ 

4.579784 

16 OC 

0.55 

9.0X10 3 

3.95424 


.575 

32.0 

4.505150 


.90 

8.0 

3.90309 


1.04 

28.0 

4.447158 






2.48 

18.0 

4.255272 


2.39 

4.0 

3.60206 


3.00 

14.0 

4.146128 


2.66 

3.5 

3.544068 


3.49 

11.5 

4.060698 


3.00 

3.0 

3.477121 






3.26 

2.76 

3.440909 

1300 

0*55 

22.0X10 3 

4.342423 

180C 

0.08 

5.0xl0 3 

3.69897 


1.16 

18.5 

4.267172 


.91 

3.55 

3.55022 


2.28 

12.5 

4.096910 


2.68 

1.65 

3.217484 


2.79 

10.0 

4.0 






3.30 

8.0 

3.903090 






3.90 

6.0 

3.778151 





1400 

0.40 

18.0*10 3 

4.255272 






.81 

14.6 

4.164353 






2.86 

7.0 

3.845098 






3.08 

6.0 

3.778151 






3.52 

5.0 

3.698970 





1500 

0.66 

ll.OxlO 3 

4.04393 






.77 

9.4 

3.973128 






1.41 

8.0 

3.903090 






2.12 

6.0 

3.778151 






2.45 

5.0 

3.698970 






3.02 

4.0 

3.602060 








TABLE III 


RUPTURE DATA FROM NIMONIC 80A 


(FROM REF. 5) 


650° 

c 

700° 

C 

750° 

c 

Stress , 

Life, hr 

Stress, 

Life, hr 

Stress, 

Life, hr 

tons/sq in. 


tons/sq in. 


tons/sq in. 


30 

274 

23 

208 

17 

138 

28 

481 

21 

443 

16 

230 

26 

898 

19 

683 

14 

419 

24 

1,292 

16 

1,735 

12 

852 

22 

2,655 

13 

4,836 

10 

1,857 

20 

5,270 

10 

10,896 

8 

4,450 

18 

8,171 

7 

34,053 

6 

13,089 

16 

13,386 



4 

22,657 




(T a , log t J 


T + 460 


T + 460 

(a) Constant-nominal- 

(a) Constant -nominal- 

(a) Constant-nominal 

stress plots. 

stress plots. 

stress plots . 


(T + 460) (log t + C) 

T - T a 

2. 3D 


log t - log t a 

te T+460 

(b) Master curve. 

■(b) Master cu've* 

(b) Master curve. 

Figure 1. - Larson- 

Figure 2 . - Li lear 

Figure 3. - Dorn param 

Miller parameter 

parameter me :hod for 

eter method for ex- 

method for extra- 

extrapolat io i (ref. 

trapolation (ref. 3) 

polation (ref. l) . 

2). 
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Figure 10. - Evaluation of best values for T a and log t a for Nimonic SOA using data from reference 








